Universality and conformal non-invariance in self-affine rough surfaces 
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We show numerically that the roughness and growth exponents of a wide range of rough surfaces, 
such as random deposition with relaxation (RDR), ballistic deposition (BD) and restricted solid- 
on-solid model (RSOS), are independent of the underlying regular (square, triangular, honeycomb) 
or random (Voronoi) lattices. In addition we show that the universality holds also at the level of 
statistical properties of the iso-height lines on different lattices. This universality is revealed by 
calculating the fractal dimension, loop correlation exponent and the length distribution exponent of 
the individual contours. We also indicate that the hyperscaling relations are valid for the iso-height 
lines of all the studied Gaussian and non-Gaussian self-affine rough surfaces. Finally using the direct 
method of Langlands et.al we show that the contour lines of the rough surfaces are not conformally 
invariant except when we have simple Gaussian free field theory with zero roughness exponent. 



I. INTRODUCTION 

Many theoretical and numerical efforts have been fo- 
cused on the study, characterization, and understanding 
of stochastic surface patterns, for various growth mod- 
els relevant to non-equilibrium processes [l|-|3|. The sur- 
face roughening phenomena have been intensively stud- 
ied via various discrete models and continuum equations. 
Scaling properties have been observed in time and space 
fluctuations of these surfaces and such interfaces show 
self-similar or self-affine properties Q. 

Various growth models are often characterized and 
classified by three exponents, the roughness exponent, 
a, the dynamical exponent, z, and the growth exponent, 
/?. Most of the work is thus devoted to identifying the 
different universality classes to which the models stud- 
ied belong [3t11|- For example various discrete growth 
models such as the ballistic deposition (BD) I9l.[l0ll , Eden 
Q and the restricted solid-on-solid (RSOS) [ll[ models 
were known to be described by the Kardar-Parisi-Zhang 
(KPZ) equation in one dimension [6[. Random depo- 
sition with surface relaxation (RDR) is another impor- 
tant discrete growth model that can be described by the 
Edwards- Wilkinson (EW) equation [J]. 

One of the most popular methods for characterization 
of the rough surfaces is based on the concept of fractal 
properties of iso-height lines called contour lines. The 
contour plot consists of closed non-intersecting loops in 
the plane that connect points of equal heights. Several 
experimental and numerical studies obtained the char- 
acterization of the fractal properties of loop ensembles 
of (2+l)Dimensional self-affine rough surfaces such as 
glassy interfaces and turbulence [l2|, (2+l)Dimensional 
fractional Brownian motion [13J, discrete scale invariant 
rough surfaces 14J], KPZ surfaces [l5j, the multi- fractal 
surfaces [la], experimental data coming from the AFM 
analysis of WO3 surfaces [17J and also STM images of 



rough metal surfaces [18j. One of the most important 
finding in this direction is the dependence of the differ- 
ent exponents i.e. the fractal dimension of one contour, 
Df, the length distribution exponent, r, and the loop 
correlation-function exponent, xi, related to the contour 
lines of mono-fractal rough surfaces to the only univer- 
sal parameter in the system, i.e. roughness (Hurst) ex- 
ponent. This conjecture confirmed by many numerical 
simulations but there is no theoretical proof yet [18| . 

Most of the discrete surface growth models in (2+l)Di- 
mension have been simulated on the square lattice with 
periodic boundary conditions with the topology of the 
torus Q. There are a few works devoted to the study 
of the dynamical scaling exponents of surface growth on 
substrates with fractal structures Qjl, [20| ■ Here we nu- 
merically study the discrete growth models including ran- 
dom deposition (RD), RDR, BD and RSOS on different 
two dimensional lattice types such as square (5), hon- 
eycomb (%), triangular )T), and also Voronoi random 
structure (V). We show that universality holds, indepen- 
dent of the lattice type. The dynamical exponents a and 
j3 for the surface growth process and also the geometri- 
cal scaling exponents Df, xi and r for the contour lines 
of such surfaces are not affected by the change of the 
substrate's lattice type. 

Another method to study two dimensional critical sys- 
tems is by Schramm-Loewner evolution (SLE K ). It gives 
a powerful tool to classify all the conformally invariant 
curves in two dimensions, for review see [21] . Recently, 
some studies argued that scaling exponents of 2D sys- 
tems i.e. zero-vorticity lines of Navier-Stokes turbulence 
[22j | and domain walls in statistical models [2l[ can be 
determined by the diffusivity constant n. We will show 
that the contour loop ensemble of RDR model, indepen- 
dent of lattice structure is conformal invariant object and 
could be described by SLE4. Our numerical calculations 
show that the contour lines of RSOS model are not con- 



formal invariant objects, however they can still possibly 
be classified by the Loewner's equation. 

The structure of this paper is organized as follows: In 
the next section we will investigate the scaling relations 
for the discrete surface growth processes. Scaling behav- 
iors and corresponding properties of contour lines are also 
presented in detail in this section. The numerical results 
to measure scaling exponents for the stochastic surfaces 
i.e. RD, RDR, BD and RSOS models on different lattice 
types are discussed in the third section. In this section we 
also identify geometrical properties of loop ensembles by 
means of conformal invariant test. In the last section we 
will summarize our studies and we draw some concluding 
remarks. 



II. SCALING OF ROUGH SURFACES 

There are many non-equilibrium surface growth pro- 
cesses which exhibit scaling properties. Different models 
with the same scaling exponents are grouped into uni- 
versality classes characterized by the same value of the 
critical exponents. Two methods to get some informa- 
tion about the scaling properties of the surface growth 
processes are i) the dynamic evolution of the aggregate 
interface and, ii) iso-height lines of the surfaces in the 
saturation regime. 



A. Dynamical scaling exponents 

The simplest quantitative behavior of a given ag- 
gregate interface is its interface width W(t, L) = 



(rdSi=i(^i(0 — M^)) 2 ) where h(t) is spatial average 
of height at time t and hi(t) is the local height variable at 
the site i. The averaged roughness over different config- 
urations show scaling behaviors. The width is saturated 
as W ex t@ with growth exponent (3 and W oc L a with 
roughness a for short time and long time limits, respec- 
tively. The scaling exponents a and /3 are used to char- 
acterize a given universality class of surface growth pro- 
cess ■ Another quantity that describes dynamic of the 
growing rough surfaces is the height-height correlation 
function C(r, t) = {[h(r + r, t) - h(r , t)]). The heights 
separated by the short distances r <C L are fully corre- 
lated and the correlation function scales as C(r, t) ex r 2a 



B. Geometrical scaling exponents 

For a given scale invariant height configuration h(r), 
at the level cut at the mean height h(r) = ho, there 
are many closed non-intersecting contour loops that con- 
nect points of equal height. A few scaling functions and 
scaling exponents need to characterize the size distribu- 
tion of such contour lines. Contour loops with length s 



and radius R are scale invariant and follows a power law 
(s) ~ (R) D} , where Df is the fractal dimension of a con- 
tour loop. The contour line properties can be described 
by the probability distribution of contour lengths P(s). 
This function measures the probability that one loop has 
length s and it follows the scaling behavior P(s) ~ s~ T 
where r is a scaling exponent. Another interesting quan- 
tity with the scaling property is the loop correlation func- 
tion G c (r). This function is a probability measure for the 
two points in the plane where separated by the distance 
|r| lie on the same contour. For iso-height lines on the 
grid with lattice constant c and in the limit |r| ^> c, 
the two point correlation function is also scale invariant 



G c (r) 



H 2 



where xi is the loop correlation-function 



exponent. For the contour loops on a Gaussian surface, 
xi = | is conjectured [241 ] . There is no any mathematical 
proof for this hypothesis yet. The estimated value xi = h 
is confirmed by various numerical works for larg e class of 
the known self-affine rough surfaces [12|, [l4|, LL8J ■ 

For the Gaussian self-affine rough surfaces with the 
roughness exponent a and the scaling exponents Df, r 
and xi, it was shown that the two well known hyper- 
scaling relations 



D f (r-1)=2 



D f (r-3) = 2xi-2, (1) 



are valid [18[. Finally there is another measure associ- 
ated with a given contour loop ensemble which is called 
the two point correlation function for contour lines with 
length s, G s . Scaling properties of the contours force 
G s (r) to scale with s and r as G s (v) ~ s m jr| _ ™/G s (-^) 
where the new exponents satisfy [12[ 



2-n= Df(l -m). 



(2) 



In the next section we numerically calculate all intro- 
duced exponents a, j3, Df, r, Xi, n and m for RD, RDR, 
BD and RSOS on different lattice structures and we will 
show the validity of the two scaling relations in Eq. ([T]) 
and also the Eq. ©. 



III. NUMERICAL RESULTS 

The lattice version of RD, RDR, BD and RSOS mod- 
els are simple to describe. In all of them particles fall 
on vertical direct line onto the substrate from a random 
position above the surface |2j. In the simple on-top site 
RD model, the particle deposits on randomly chosen site 
position. In RDR model a particle is released on the top 
of a selected column, but it docs not stick to the surface 
immediately and it can diffuse to a nearest neighbor site 
of lower height [2, |j| ■ Surface diffusion generates corre- 
lation between neighbor sites and it causes surface width 
saturation effect. In BD model, falling particle sticks to 
the old one either on the top or a nearest neighbor oc- 
cupied site [2, Hi [111 ■ In RSOS surface growth process 
the particle aggregates to the substrate if and only if the 
difference between heights of all pairs of nearest-neighbor 



columns satisfies the condition Ah < 1; otherwise, if this 
condition is not met, the corresponding aggregation at- 
tempt is rejected from the system fill ]. 

In order to find dynamical and geometrical scaling 
properties of RD, RDR, BD and RSOS models on the 
initially flat substrates, we have generated these mod- 
els on different lattice types, i.e. honeycomb, triangular, 
square and Voronoi with size L = 512 and all measure- 
ments are made using an ensemble of 2000 realizations 
(for BD model the simulation were also done on lattices 
with size L = 2048 with 500 samples). In order to reduce 
the errors due to the substrate's boundaries, we used pe- 
riodic boundary condition during particle deposition in 
horizontal and vertical directions. Each time step is de- 
fined as the number of particles to fall up surface on the 
average, which is equal to L x L. To extract the contour 
lines of the saturated surfaces at the mean height, ho we 
used Hoshen-Kopelman algorithm [25j . 



A. Distribution of heights 

To check that the surface that we have is a Gaussian 
surface or not we used the definition of local curvature 
at r and at scale b which is 12411 



M 



C b (r) = Y, IH? + be m ) - h(r)} 



(3) 



where the sum of e m 's arc a fixed set of vectors summing 
to zero. The distribution of the local curvature in Gaus- 
sian surfaces is Gaussian. In Fig. 1 one can see that the 
distribution is Gaussian for RDR but non- Gaussian for 
RDR and BD. To quantify this claim we depicted in Fig. 



2 the third 



(gkWD 



and fourth moments 



(gb(rD 



of the 



<C b (r)2)3/2 —™™ <C b (r)2)2 

local curvature which should be and 3 for a Gaussian 
surfaces respectively. Based on these two figures we con- 
clude that among the surfaces that we are going to study 
just the RDR is a Gaussian surface. It is worth men- 
tioning that by increasing the scale b, the results of RDR 
and BD surface growth models converge to the Gaussian 
process. We also repeated these computations for the 
other lattice types and the results are the same as that 
of for square lattice type. Finally we also checked the 
self- affinity of our surfaces by showing that for all the 
surfaces we have ([Cf,(r)] 9 ) ~ b qa . 



B. Scaling exponents 

One approach to show that two rough surface mod- 
els belong to the same universality class is to compute 
all scaling exponents of the models. In this sub-section 
we report measured values of the scaling exponent for 
various models. 

a. The growth scaling exponent 
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FIG. 1: (Color online) The distribution of the local curvature 
in RDR, BD and RSOS model. To make more sense we shifted 
various plots vertically. 



In Table. HI we have shown the results for various 
growth models on different lattice types. Our measure- 
ments of the scaling exponent j3 on the square lattice are 
consistent with previous studies @, [26( . 

Our result confirms logarithmic scaling of roughness 
W(t, L) ~ ln(t) for RDR model where is in agreement 
with the predictions [2], |4| . 

b. Roughness exponent 

To measure the roughness exponent a we have cal- 
culated the height-height correlation function C(r) with 
respect to the distance separation r for saturated time. 
Using numerical calculations for small values of r, the 
roughness exponent can be read where Table [TT| reports 
all the measured values of the roughness exponent for 
BD and RSOS lattice growth models on different lat- 
tice types. Results from numerical computation of a for 
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FIG. 2: (Color online) The third and fourth moments of the 
local curvature in RDR, BD and RSOS models on square 
lattice type. 



TABLE I: Numerical values of the scaling exponents /? for 
different surface growth models on the different lattices (H: 
Honeycomb, <S: Square, T: Triangular, V:Voronoi). The 
numbers inside the parenthesis are the error bars of the last 
digits. 



Lattice type 



Model H 



T 



V 



RD 0.500(1) 0.500(1) 0.500(1) 0.500(1) 
RSOS 0.240(2) 0.240(5) 0.240(5) 0.240(3) 
BD 0.24(1) 0.24(1) 0.23(2) 0.24(1) 



TABLE II: Numerical values of the scaling exponents a for 
BD and RSOS models on different lattice types. 



Lattice type 


Model U 


S T 


V 


RSOS 0.400(2) 
BD 0.30(4) 


0.400(2) 0.400(3) 
0.29(4) 0.29(4) 


0.400(5) 
0.30(4) 



TABLE III: Numerical values of different geometrical expo- 
nents xi , Df and t — 1 at la confidence interval for the discrete 
surface growth models on different lattices. 









Lattice type 






H 


S 


r 


V 




RDR 


1.0(1) 


1.0(1) 


1.0(1) 


1.0(1) 


2xi 


RSOS 


1.0(1) 


1.0(1) 


1.0(1) 


1.0(1) 




RD 


0.53(5) 


0.47(5) 


0.49(4) 


0.52(5) 




BD 


1.0(2) 


1.0(2) 


1.0(2) 


1.0(2) 




RDR 


1.50(1) 


1.50(1) 


1.49(1) 


1.49(2) 


D f 


RSOS 


1.30(1) 


1.30(1) 


1.30(1) 


1.30(1) 




RD 


1.73(1) 


1.73(1) 


1.75(1) 


1.73(1) 




BD 


1.36(2) 


1.36(2) 


1.36(2) 


1.36(2) 




RDR 


1.33(1) 


1.34(1) 


1.34(1) 


1.33(1) 


T-1 


RSOS 


1.22(1) 


1.23(1) 


1.22(1) 


1.22(1) 




RD 


1.18(1) 


1.19(1) 


1.19(1) 


1.19(1) 




BD 


1.32(2) 


1.31(2) 


1.32(3) 


1.32(2) 




RDR 


0.49(2) 


0.50(2) 


0.50(2) 


0.51(2) 


n 


RSOS 0.69(2) 


0.69(2) 


0.68(2) 


0.69(2) 




RD 


0.30(2) 


0.29(2) 


0.30(2) 


0.31(2) 




BD 


0.66(2) 


0.67(3) 


0.65(2) 


0.66(2) 




RDR 


0.02(3) 


0.02(2) 


0.02(2) 


0.02(2) 


m 


RSOS 0.03(3) 


0.03(3) 


0.03(3) 


0.02(3) 




RD 


0.03(3) 


0.02(3) 


0.02(2) 


0.03(3) 




BD 


0.04(5) 


0.04(4) 


0.03(4) 


0.03(3) 



growth processes on the square lattice is in good agree- 
ment with the previously reported values @, [26[ . 

For RDR model deposited on different lattice types, 
we observed C(r) ~ ln(r) [2,|j]. 

c. Geometrical scaling exponents 

The most important exponent in fractal contour lines 
is the loop correlation function exponent xi . For a given 
loop ensemble we followed the algorithm reported in the 
Ref. [12| to find the correlation function G c (r). The first 
part of the Table IIIII shows that the measured values for 
the exponent Xi for discrete growth processes i.e. RDR 



model with the Gaussian statistics are in agreement with 
the predicted value 2xi = 1, within the statistical error. 
According to our results, it is quite interesting that the 
relation 2.ti = 1 is valid also for non-Gaussian interfaces 
such as RSOS model. We observed 2x; = i for the RD 
model which is in agreement with the predicted value 
for the two point correlation exponent of the percola 



tion model [27j. The difference between RD model and 



the other growth processes comes from the correlation 
between lattice sites in the deposition process. There 
are no any correlations between nearest neighbor sites 
in the RD model. The same calculations are done for 
the BD model on different lattice types. The large er- 
ror bars which lead to the less agreement in the equality 
2xi = 1 come from finite size effects. One can com- 



pute the perimeter, s, and gyration, R, radius of con- 
tour loops to calculate fractal dimension Df. Here R is 
defined by i? 2 = ± ^ =1 [fa - x L f + ( Vi - y c ) 2 ] , with 

x c = ~N X)»=i x i and Vc — 7f J2i=i Hi- In tne second part 
of the Table UTTl we have shown the measured values of the 
fractal dimension of contour loops for all surface growth 
processes on different lattice types. The next remark con- 
cerns the probability distribution of contour length P(s) 
with the scaling exponent t. We measured r for differ- 
ent models, results are shown in the third part of Table 
IIII1 A suitable value of the parameter r obtained form 
the best fit to ln(P(,s)) versus ln(s) (We measured r — 1 
instead of r). Finally we also calculated G s (r) and the 
exponents m and n; the results are listed in Table 1IIII 

d. Scaling relations 

Finally the relation between measured scaling expo- 
nents a, Xi, Df, r, 77i and n given by Eq. [T]and Eq. [5] 
can be examined. We have depicted the results for all 
growth models on different lattice types in the Table IIVI 
As can be seen from the Tables. U [XT] and IIIII the error 
bars for the BD results are very large. It seems that to get 
conclusive results one needs to simulate the BD model on 
very large lattice sizes j28| which is very time consuming. 
However, due to the finite size effect, we concluded that 
our results are roughly consistent with the predictions. 



C. Conformal invariance (CI) test 

It is well- known, for a recent discussion see [29|, that 
the action corresponding to Gaussian self-affine rough 
surfaces is conformally invariant just for the surfaces with 
zero roughness exponent. However, recently many au- 
thors claimed that despite the non-conformal invariant 
height ensembles of these surfaces, their contour lines 
might be conformally invariant 15, 26, 30]. In the follow- 
ing we examine this guess by using the direct conformal 
invariance test used first in |3l| . Suppose a conformal 
mapping defined by w = gt(z). For any domain T> with 
boundary C, one can find a conformal map g t where maps 
T> to V and C to C. For the critical statistical systems, 
the measure on distributions on C is obtained by trans- 
port of the measure on distributions on C using conformal 
mapping g t from T> to T>' where this measure is invariant 
at the critical point of the model. 

Statistical systems such as, percolation and Potts mod- 
els at the critical point, exhibit a unique spanning cluster, 
where connect the boundaries of domain from one side 
to the other (i.e. left to right or up to down) along the 
spanning cluster |32j. In the scaling limit when the size 
of the system goes to infinity, one can define the crossing 
probability nh (7r„), the probability of a system to per- 
colate only in the horizontal (vertical) direction [31] . In 
the conformal invariant statistical systems, the measure 
i7i (jTv) is unchanged under any conformal map gt where 
it maps a given boundary C to C J3l|, I33W35J . 



The conformally invariant curves should be statisti- 
cally equivalent to Schramm Loewner evolution (SLE K ). 
In this method the conformal map gt{z) which maps 
the half-plane minus the trace "ft into itself, obeys the 
Loewner's formula dtgt = 2/(gt — £t) where the driving 
function £ t is related to the Brownian motion with the 
diffusivity k. The fractal dimension of such scaling curves 
is given by the relation df = 1 + k/S. In [15, 26J,[30J, using 
a discrctized Loewner equation and iterative conformal 
slit map they extracted the drift £4 for iso-height lines in 
the different rough surfaces. They found that an ensem- 
ble of the driving function £ t for these models looks like 
converging to a Gaussian process with zero mean (£ t ) w 
and the variance (£ 2 ) « Kt with k compatible with the 
fractal dimension of the curves. It shows for example that 
the model with logarithmic height correlation can be re- 
lated to SLE4. However, it seems that the calculations 
based on Loewner evolution (using Schramm Loewner 
equation to find k) are not accurate enough. Especially 
it is not easy to find the exact Gaussian distribution for 
the increments of the driving function and also it is dif- 
ficult to show stationarity of increments of the drift to 
conclude that the drift is really a Brownian motion. For 
example even for RDR (which is a Gaussian free field) we 
know that usually the drift is not just a Brownian mo- 
tion, it is a complicated function related to the Brownian 
motion [36|. 

An alternative way is to apply directly the confor- 
mal invariant test to a two dimensional slice of ran- 
dom rough surface. For a given rough surface sam- 
ple h(x,y), a horizontal cut is made at a certain level 
h c = (h) +5y/([h — (h)] 2 ), where the symbol (...) is aver- 
aging over ensemble of rough surfaces. A set of nearest- 
neighbor connected sites of positive (negative) height 
form a cluster. The cluster is spanning one if spans two 
opposite sides of a given domain, such that infinite con- 
nectivity first occurs, and one can measure the crossing 
probabilities irh{8) and tt v (S) for a given domain geome- 
try. 

In this study we considered square boundary as C and 
rhombus domains with different angle as C . There is 
a conformal map that takes interior of our choice of C 
to the interior of C' and takes vertices to vertices and 
sides to sides. The basic method to perform this test 
composed of three steps: (I) draw the boundary C ( C) 
defining cluster ensemble on the lattice. (II) assign a 
state +1 (—1) to each site of the lattice located inside C 
(C) when h(x 7 y) > h c (h(x,y) < h c ). (Ill) find cluster 
of nearest-neighbor connected sites with the same values 
(for example +1) by using Hoshenn-Kopelman (HK) al- 
gorithm [23] ■ These three steps are repeated for all of the 
configurations of level set ensemble for various <j's to find 
the largest cluster where spans left and right, to measure 
iTh (up and down to measure n v ). The expected value of 
7T/t (tTii)) is then the ratio of the number of configurations 
spanning two opposite horizontal or vertical sides of the 
boundary C ( C) to the sample size. 

We examined the above test for two dimensional slice 



TABLE IV: Verification of three basic hyper-scaling relations for discrete surface growth process on different lattice types. 







D,(t-1) 


2-a 


Df(r- 


3) 2-2xi 


D f (l-m) 1-n 




H 


2.04(3) 


2 


1.42(3) 


1.47(5) 


1.68(4) 


1.70(2) 


RD 


S 


2.06(5) 


2 


1.40(5) 


1.53(5) 


1.69(4) 


1.71(2) 




T 


2.08(6) 


2 


1.42(6) 


1.51(4) 


1.71(3) 


1.70(2) 




V 


2.06(5) 


2 


1.40(5) 


1.48(5) 


1.68(4) 


1.69(2) 




H 


1.79(6) 


1.70(4) 


0.93(4) 


1.0(2) 


1.30(6) 


1.34(2) 


BD 


S 


1.78(8) 


1.71(4) 


0.94(6) 


1.0(2) 


1.30(6) 


1.33(3) 




r 


1.79(6) 


1.71(4) 


0.93(4) 


1.0(2) 


1.32(6) 


1.35(2) 




V 


1.79(6) 


1.70(4) 


0.93(4) 


1.0(2) 


1.32(5) 


1.34(2) 




H 


2.00(3) 


2 


1.01(3) 


1.0(1) 


1.47(4) 


1.51(2) 


RDR 


s 


2.01(3) 


2 


0.99(3) 


1.0(1) 


1.47(3) 


1.50(2) 




T 


2.00(3) 


2 


0.98(3) 


1.0(1) 


1.46(3) 


1.50(2) 




V 


1.98(4) 


2 


1.00(4) 


1.0(1) 


1.46(3) 


1.49(2) 




H 


1.59(2) 


1.60(1) 


1.01(2) 


1.0(1) 


1.26(4) 


1.31(2) 


RSOS 


S 


1.60(2) 


1.60(1) 


1.00(2) 


1.0(1) 


1.26(4) 


1.31(2) 




r 


1.59(2) 


1.60(1) 


1.01(2) 


1.0(1) 


1.26(4) 


1.32(2) 




V 


1.59(2) 


1.60(1) 


1.01(2) 


1.0(1) 


1.27(4) 


1.31(2) 



of RD, RDR, RSOS and also synthetic (2+l)Dimensional 
self-affine surfaces which known as fractional Brownian 
motion (2D FBM) [13[ for different 8's and we measured 
crossing probability iTh(6) (n v (8)) for different boundaries 
C (square with side length L) and C (rhombus with angle 
9 and side L). In order to have an ensemble of surface 
height profiles, we have obtained 2 x 10 5 samples from the 
grown rough surfaces with side L = 512 for each model. 
We measured the probabilities irh(S) and x v (8) that at 
each level height 8, an infinite island spans two opposite 
boundaries of the domains. 

As shown in Fig [3] the measured crossing probabilities 
for RDR model on rhombus domains with different angles 
6, cross each other at 8 = S c , implying that the crossing 
probabilities remain constant by conformal transforma- 
tion. This means that the two dimensional slices of RDR 
model at the critical level 8 C is possibly conformally in- 
variant and the crossing probabilities do not change by 
conformal transformations. We also applied this test to 
2D FBM with Hurst exponent < H < 1 and we ob- 
served conformal invariance only for H — 0. We observed 
that the 2D FBM with H > is conformal non-invariant 
and the spanning probability irh (8) do not show any fixed 
point. In addition we presented in Fig. [3] the curves 
7T/i(<5) obtained for RSOS model on different boundaries. 
It shows that the crossing probability for two dimensional 
slice of RSOS model for different 8's, change by conformal 
transformation of the boundary. Increasing the number 
of ensembles and the size of the system did not make sig- 
nificant change in this pattern. We concluded that this 
model is also conformally no n- invariant. 

Our results show that the contour ensemble of RD 
model on the triangular lattice at the mean height (8 = 0) 
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FIG. 3: (Color online) The spanning probability TVh{5) for 
two dimensional slice of RSOS model on rhombus domain 
with different angles 6. Inset: The same quantity for RDR. 



is conformally invariant. Two dimensional slice of RD 
model at the mean level for other lattice types is con- 
formal non-invariant. This is obvious if we think to the 
clusters of RD as a simple percolation clusters. Our nu- 
merical results are summarized in Tabel. fV)) . 



IV. CONCLUSION 

In this paper we showed that many properties of 
rough surfaces such as the roughness exponent and 
dynamical exponent are independent of the underlying 
lattices, see Table. Q] and [TTJ We also showed that the 
fractal dimension and correlation exponent and many 
other exponents of contour lines of the surfaces are also 



TABLE V: CI test results of two dimensional loop ensembles 
of random rough surfaces. 





Lattice type 


model 


S T 


RD 


No Yes 


RDR 


Yes Yes 


RSOS 


No No 


FBM(H ^ 0) 


No No 



not conclusive we were able to show that the extracted 
critical exponents of the contour lines of RDR, RSOS 
and RD follow the hypcrscaling relations independent of 
having a Gaussian or non- Gaussian surface, see Table. 
IVl Finally using direct conformal maps we showed that 
the contour lines of self-afnne rough surfaces (RSOS, BD 
and fractional Gaussian rough surfaces) with non-zero 
roughness exponents are not conformally invariant. To 
be specific just the contour lines of RDR on different 
lattices and RD on the triangular lattices are conformally 
invariant, see Table. IVl 



independent of the underlying lattices and the distribu- 
tion of the heights. All the calculations were carried out 
on four different regular and irregular lattices such as 
the square, triangular, honeycomb and Voronoi lattices, 
see Table. IIIII Although the calculations for BD were 
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